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Abstract 

The electron transport in a ID conductor with an isolated local defect such as an impurity or a 
non-adiabatic contact is studied theoretically. New regime of conduction in correlated ID systems is 
predicted beyond the well-known regime of tunneling resulting in the power-law I-V-curves. In this 
regime a quantum wire becomes "opened" at voltage biases above the threshold value Vr, giving rise 
to a rapid increase of the dc current, /, accompanied by ac oscillations of frequency / = I/e. The 

effect of ac generation is related to sliding of the Friedel oscillations of electronic density produced 

O 

by the defect. Manifestations of the effect resemble the Coulomb blockade and the Josephson 

effect. The spin bias applied to the system is shown to affect the I-V curves due to violation of 

^y-s the spin-charge separation at the defect site. In short quantum wires of length L < Lq ~ hvp/eVr 

{\j and at high temperatures T > Tq ~ eVr/fce the effect is destroyed by fluctuations. The threshold 

in 

_J voltage Vr is determined by 2/cir-component of impurity potential renormalized by fluctuations. 

The ID conductor is described in terms of the Tomonaga-Luttinger Hamiltonian with short range 
• • or long-range Coulomb interaction by means of the bosonization technique. We derive boundary 

conditions that take into account relaxation in the leads and permit to solve non-equilibrium 
problems. Charge fluctuations are studied by means of Gaussian model which can be justified 
strictly in the limit of large voltages or strong enough inter-electronic repulsion. Spin fluctuations 
are taken into account strictly by means of the refermionization technique applicable in case of 
spin-rotation invariant interaction. 

PACS numbers: 73.63.Nm,73.21.Hb,73.63.-b, 72.15.Nj, 



I. INTRODUCTION 

Basic electronic properties of three-dimensional (3D) solids are usually well described 
within the Landau's Fermi-liquid picture where low-energy excitations are single-electron 
quasiparticles that in many respects behave like non-interacting electrons. In the Fermi- 
liquid the interaction modifies energy spectrum of electrons and makes their life-time finite 
but still very long if the energy of electron is near the Fermi energy. This is not the case 
in ID systems where interaction in some sense is always strong and cannot be considered 
as a small perturbation. The picture of single-electronic excitations fails in correlated ID 
systems, and the only low energy excitations are bosonic charge and spin collective modes 
with a sound-wave-like spectrum. These excitations form the Luttinger liquid (LL) that is 
an alternative to the Fermi liquid for ID electronic systems (for a review see Ref. Q] and 
EJ). Such distinctive features of the LL as spin-charge separation and power-law suppression 
of tunneling into ID systems have been confirmed experimentally, see e. g. Ref. El There 
are different realizations of ID electronic systems demonstrating properties of the LL. The 
examples are semiconductor-based quantum wires in which dimensionality of the conduction 
electrons is reduced by dimensional quantization^, metallic linear chains on Si surfaces^, 
carbon nanotubea^, conducting polymers^, and edge states in FQHE, and 2D topological 
insulators 9 . There are also evidences of the LL-like effects in electronic transport in strongly 
anisotropic quasi-lD conductor ^ 10 1 11 1 where the LL state can be stabilized by defects^^ or 
by formation of the charge-density wave (CDW) gap induced by electron-phonon coupling^. 

Inter-electronic interaction greatly affects electronic transport in conducting ID systems. 
In particular, the back-scattering component of the impurity potential in ID systems with 
repulsive inter-electronic interaction (interaction parameter K p < 1) scales to infinity under 
renormalization group transformations. Hence, even isolated impurities form effectively large 
barriers and strongly suppress conductance^"^. This is reflected in power-law dependencies 
of conductance on voltage and/or temperature. This was confirmed experimentally in various 
ID systems, including semiconductor quantum wires^ and carbon nanotubes™. 

Strong effect of a single impurity on the electronic transport in ID electron system is 
related to Friedel oscillations induced by impurity^, and the power-law I-V curves were de- 
scribed in terms of macroscopic tunneling between different minima of a washboard potential 
tilted by external voltage bias. The periodic potential is associated with interaction between 



smooth electronic density of the wire and Friedel oscillations. On the other hand, the limit 
of strong interaction between electrons in solids usually leads to the Wigner crystallization. 
But it is well-known that purely one-dimensional (ID) crystals cannot existP^tLSJ because the 
long-range order is destroyed by fluctuations. So, strictly speaking, ID Wigner crystals do 
not exist neither. However, in ID gas with Coulomb repulsion density-density correlation 
functions contain oscillatory parts: the first one is at wave vector 2/c^p and decays approxi- 
mately as 1/x, and the second one is at wave vector Akp and decays extremely slowly 2 ^, like 
g-cvinx^ ^at is slower than any power-law. As the period corresponding to Akp oscillations 
is exactly the average inter-electron spacing, such a system can be considered as a ID Wigner 
crystal with pseudo-long-range order^. In case of short range inter-electronic interaction 
(which takes place in gated quantum wires where the long-range part of the Coulomb inter- 
action is screened by electrons in the metallic gate) the density correlations both at 2kp and 
at 4kp decay as power-law, with the exponent of 4kp correlations being smaller. Despite 
the fact that there are oscillating terms in the correlation function, expectation values of 
these terms are absent in homogeneus LL as they are washed out by fluctuations. But the 
oscillating terms are present near defects since they strongly pin the LL and suppress density 
fuctuations. Sliding of the electronic crystals is known to contribute to conductivity of the 
system, the most studied case being quasi-lD CDW compounds 21 . As defects violate the 
space invariance they pin the CDW and for the CDW sliding the driving electric field must 
be larger than a threshold value. Above the threshold field the CDW starts to slide and to 
contribute to the current resulting in non-linear collective electron transport and generation 
of the so-called narrow-band noise at washboard frequencies^ corresponding to a shift of the 
CDW by one period. As long as the LL can be interpreted as a ID form of the ID Wigner 
crystal, one can expect that the similar dynamic regime of depinning and sliding must ap- 
pear in correlated ID electron system as well. Thus, one can expect that at large enough 
voltages the system can roll out from the minimum of the washboard potential, namely, 
at voltages exceeding a threshold value determined by the slope of the washboard. Such a 
behavior is expected, at least, in the quasiclassical limit when quantum fluctuations at the 
impurity site are suppressed by strong interaction between electrons. This corresponds to 
very small values of the interaction parameter K p in case of short-range interaction or to 
Coulomb interaction when effective K p depends on momentum and is strongly suppressed 
in the long wave-length limit because of the long-range nature of the Coulomb interaction^. 



Such a scenario was addressed in our letter^ where the dynamic regime of conduction in 
the LL was predicted for quasiclassical case of K p <C 1, in this regime the dc current is 
accompanied by oscillations with the washboard frequency / = I/e. 

Full I-V curves of a single- channel LL with a single impurity were studied by means of 
thermodynamic Bethe ansatz technique by Fendley et aP. Egger and GrabertP^l calculated 
the I-V curves for specific value of interaction parameter K p = 1/2 using the refermioniza- 
tion technique which makes the Hamiltonian quadratic and, hence, solvable exactly. Cal- 
culation^ of the current in the limit of strong interaction, K p <C 1, resulted in S-shaped 
I-V curves, and the authors concluded that by putting such a quantum wire into a properly 
designed load circuit, self-sustained current oscillations can be generated. These oscillations 
are a consequence of instability induced by negative differential conductivity, and their origin 
different from the mechanism discussed in Ref. 22 and considered here. 

In all these papers the equilibrium distribution function of incident particles (non- 
interacting fermions, kinks and antikinks, etc) was assumed while distribution of the trans- 
mitted particles was not the equilibrium one. However, bosonic excitations of the LL are 
reflected from the leads even in case of adiabatic contacts because the interaction is present 
inside the wire and is absent in the leads. Reflection coefficient equals r = ~ p which 
approaches to unity in the limit K p — > 0. As long as relaxation processes inside the wire 
are not taken into account in the Hamiltonian, the transmitted waves do not relax inside 
the wire even in the limit of very large length. Thus the incident waves must consist in part 
of the waves reflected from the contact. So distribution function of incident particles must 
not be necessarily the equilibrium one, and this applies equally to fermions derived from 
bosons after the refermionization. In this case one needs to take into account relaxation 
processes induced by coupling of the quantum wire to 2D or 3D Fermi liquid of the leads 
which can be considered as a heat bath. This can be made by means of boundary condi- 
tions. They allow to find non-equlibrium fluctuations inside the wire in terms of equilibrium 
density fluctuations in the leads. One of our goals is to include such effects into study of 
the electronic transport in quantum wires. This will be performed using boundary con- 
ditions for operators which generalize the radiative boundary conditions derived by Egger 
and Grabert 24 ' for expectation values. Furthermore, using boundary conditions for non-ideal 
contacts, we demonstrate that Friedel oscillations which appear near non-adiabatic contact 
result in non- stationary effects similar to those of impurities. Note that we consider non- 



ideal contacts in which backscattering is induced by the geometry, and this is different from 
the model of Ref. [2H1 in which backscattering was derived in a microscopic lattice model 
with a spatially dependent nearest-neighbor interaction turned on near the contact from 
zero to a bulk value. In the local LL model characterized by a spatial dependence of the LL 
parameters the backscattering does not appear^"^. 

The dynamic regime of conduction was considered in Ref. 22 in the limit of strong inter- 
electronic repulsion, when fluctuations at the impurity site are small. The results were also 
limited by the case of very long conducting channels, longer than the damping length of ex- 
citations due to coupling of electrons inside the wire to a dissipative bosonic bath (phonons, 
density fluctuations in a metallic gate, and so on) when one can neglect that current pulses 
generated at the impurity reflect from the contacts. Here we consider conduction in more 
realistic case of shorter conducting channel when fluctuations in the channel are governed 
by boundary conditions. We study solutions at various values of the applied voltage con- 
centrating on the limit of large voltages/currents, where Gaussian fluctuations dominate, 
and non-Gaussian part of fluctuations is smaller. This enables us to calculate I-V curves, 
the amplitude of the ac current and noise spectrum more strictly. The I-V curves that we 
found differ from the results of Ref. [23] and [2H In addition, we consider both a single mode 
(spin-polarized) and a spinful ID electronic systems. The interaction between electrons is 
assumed to be strong enough and of two types, either a short range characterized by a con- 
stant interaction parameter K p < 1, or the long range Coulomb interaction. The short-range 
interaction describes gated quantum wires where the long-range part of the interaction is 
screened by electrons of 3D gate electrodes. 

The structure of the paper is as follows. In Sec. [IT] we formulate the problem, and derive 
boundary conditions for equations of motion for Heisenberg operators of the displacement 
fields of the LL. Using the boundary conditions we derive equations of motion for the dis- 
placement field at the impurity site. These equations resemble equation of motion of a 
quantum pendulum or equation for the superconducting phase in the Josephson junctions. 



In Sec. [ITT] we use our equations to study electronic transport in spinless LL. Using the 
Gaussian model to account for fluctuations, we calculate I-V curves in several limits, an- 
alyze noise spectrum, study non-Gaussian corrections, and finally discuss validity of the 
Gaussian approximation. In Sec. |IV| we consider the spinful LL with strong enough inter- 
action between electrons (either very strong short-range or moderate long-range Coulomb 



interactions). This is the limit in which charge fluctuations at the defect are small. Though 
spin fluctuations are still large, they are taken into account strictly by means of refermioniza- 
tion method in spin sector. This method works in case of spin-rotation invariant interaction 
{K a = 1, interaction does not depend on spin). In Sec. Mwe show that non-ideal contacts 



induce non-stationary effects similar to those induced by impurities. In Sec. |VI| we formulate 
conclusions. 

Below we set e, h and fc# to unity, restoring dimensional units in final expressions when 
necessary. 

II. GENERAL FORMULATION 

A. Problem formulation 

We consider a correlated ID conductor containing a defect located at x — 0, and adia- 
batically connected to ideal Fermi-liquid reservoirs at x = ±L/2, or an ideal ID conductor 
without impurity connected to reservoirs by non-adiabatic contacts. The Hamiltonian of 
the system with impurity consists of two terms 

H = H + Hi. 

The first one is the bosonised Tomonaga-Luttinger (TL) Hamiltonian that maps the ID 
system of interacting electrons to free massless bosons described in terms of the displace- 
ment fields $ p (t, x) and <t a (t, x) and the conjugated fields Il„(£, x) = d x O u /n describing the 
momentum density, [& u (t,x), Il v i(t,y)] = i8(x — y)S u y. Here v = p,a denotes charge and 
spin channels, correspondingly. These fields are related to fermionic electron operators for 
right (r = +1) and left (r = —1) moving electrons with spin s 

i f A* it 1 nr 

^ rs (x) = lim e ^ 7)rse -^[© P W-^W+«(^W-^W)] ; (1) 

where i) r ^ s are ladder operators, and a^kp 1 is the small cut-off length of the Tomonaga- 
Luttinger model. The standard TL Hamiltonian in the Fourier transformed form reads^l 

Here the LL parameters K u , playing the role of the stiffness coefficients of the elastic string 
described by Hamiltonian (J2J), are related to the electron-electron interaction potential, and 



measure the strength of interaction between electrons. In the spin-rotation invariant case 
considered in our study, K a = 1, and 

K ^ = JTTm- (3) 

V 7™F 

where g(q) is the Fourier transformed interaction potential. In case of infinite ID gas with 
long-range Coulomb interaction described by the approximate form Vc(x) = , % , 2 , where 
e is a background dielectric constant and d is a diameter of quantum wire, one obtains 
g(q) = 2^K (M|)]. Thus, 

M 1 2e 2 2 /_^\ 1 

P[q) ^/l+ 7 Ko(M)' 7 ^jpc 1377T UJ ? l j 

where 7 is dimensionless parameter which measures the strength of the Coulomb repulsion 
between the electrons. 

This case was described by Schulz^ as a "ID Wigner crystal" because the system exhibits 
a 4fcp electron density modulation characterized by extremely slow (slower than any power 
law) decay of the density correlations, which is induced by long-range nature of Coulomb 
interaction. 

In case of the short-range interaction the dependence of g on wave-vector q is usually 
neglected, then K p < 1 corresponds to repulsive and K p > 1 to attractive interaction. 
The model with short-range interaction is applicable to gated quantum wires where the 
long-range part of interaction is screened by 3D gate electrodes. For quantum wires one 
can roughly estimate the value of K p within a simple model of the nanowire of diameter d 
located at distance d\ from the metallic gate 

„_ 2 1 2e 2 d + d 1 1 / c \ 2 d + d x 

K / « H In — = H — — In — . 

p TrnvF^ d 137 \vfJ vre d 

Eigenmodes of free electronic system described by Eq. [2] are bosonic waves with a sound- 
like spectrum, cu u (q) = qv u (q) where v v (q) = Vp/K u (q) is the velocity of charge [y = p) and 
spin [y = a) excitations. 

In case of the long-range interaction and finite length of the system located at \x\ < L/2 
we have to modify the Coulomb potential in order to take into account screening of the 
interaction potential by metallic leads. The exact form of the screening depends on the 
geometry of the system. We will consider 3D metallic leads forming sheets of a plane 
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capacitor connected by the quantum wire. Then the screening by the leads can be depicted 
in terms of the image charges, and the interaction potential between charges located at x 
and x' is described as 

oo 

V(x, x') = Yl [ y c{x - x' + 2nL) - V c {x + x' + 2nL + L)] , (5) 

n=— oo 

where the term with n = describes the direct Coulomb interaction, and other terms are 
induced by image charges. Its contribution to the v = p term in the Hamiltonian (J2J) in the 
coordinate representation reads 

dxdx {d x $ p (x)V(x, x')d x ^ p (x)} (6) 

Since the operator of the particle density is given by expression p = — (\/2/7i)d x <b(x), this 
term has rather transparent physical meaning. 

Interaction with the impurity is described in terms of the phases field $„(£,x) at the 
impurity position x = CP^ 

r 1 
Hi = - dx- W5{x) cos V2<S> P cos V2<t> a , (7) 

where the impurity strength W is related to the back-scattering part of the impurity potential 

Ui{x) 

W= [— e 2ikFX Ui{x). 
J a 

The forward scattering is not included because it can be eliminated from the problem by 
redefinition of the field $ p ^. Note that the form (JTI) of the impurity Hamiltonian follows 
from bosonization relations (fTl) and is related to 2 kp- components of electron density. The 
Luttinger model does not contain higher harmonics, but in more general model the higher 
harmonics are presentP. For example, for a model with non-linear electron dispersion, more 
general bosonization relations are usecP^, and the higher harmonics do appear. However, 
such harmonics do not result in qualitative effects in electronic transport studied here. 

Current and charge density perturbations in the system can be calculated in terms of & p 
by means of thermodynamic averaging of the expressions for corresponding operators 

I=^-dt%, p=-^d x %. (8) 

71 71 

Further, the displacement field determines also the operator of 2/c^-component of charge 
density perturbations 

p2k F = — cos(2/c F x + V2$ p ) cos V2$ a . (9) 
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In pure infinite ID electronic system fluctuations of the phase field $ are infinite, ($ 2 ) = 
oo, therefore, the expectation value of the right-hand side of Eq. ^ and, hence, 2kp- 
component of charge density is zero. However, fluctuations of $ can be finite at the impurity 
position, and this results in Friedel oscillations, i. e., in Ikp modulation of charge density, 
P2k F — (f>2k F )- The amplitude of oscillations decays with the distance. The larger electronic 
repulsion is, the slower Friedel oscillations decay. According to Eq. Q8J) , the current flow 
is related to time variations of $, therefore, according to Eq. (N) the current passing the 
impurity implies a space shift of the Friedel oscillations by the current. The drag of the 
Friedel oscillations by current illustrates the origin of strong suppression of conductance by 
impurities and its dependence on inter-electronic interaction. 

According to Eq. rt8J) the macroscopic current in the system can be expressed in terms of 
the expectation value of the displacement field. The latter can be found from equation of 
motion for the Heisenberg operator & p (t, x). Commuting $ p with the Hamiltonian for the 
case of short range interaction we find 

(vpl - df) $ p {t,x) = V2nv F W smV2$ p (t,x) cosV2$ (7 (t,x)5(x), (10) 

where v p (q) = vp/K p (q) is, again, the velocity of charge (plasmonic) excitations. Equation 



of motion for the spin field has similar form, it can be obtained from (10) by substitution 
subscripts p by a and vice versa. 

At contacts we apply the boundary conditions for the operator $, which takes into account 
injection of electrons induced by external bias and relaxation processes induced by coupling 
of the quantum wire to 2D or 3D Fermi liquid in the current leads. The issue on the boundary 
conditions is considered in details in the next subsection. 

B. Boundary conditions 

Boundary conditions for the single mode (spinless or spin-polarized) wire contacting with 
a 2D or 3D leads were derived in Ref. E31 Here we generalize this result for the spinful case. 
In order to derive the boundary conditions we use the ideas of the scattering approach (for 
a review see Ref. I3lj) and consider the region of the wire at |x| < L/2 as a scatterer. 

Without the loss of generality we assume that longitudinal (along the x-axis) and trans- 
verse motions are separable. The longitudinal motion in the leads is characterized by wave 



vectors k, spin s and energy Ei = j-. The transverse motion is described by energy e n , 
the total energy being e = ei + e n , where n is an index labeling transverse energies. We 
assume that electrons in the leads do not interact. Then we solve an equation of motion for 
electronic field operators in the leads using the continuity of both the field operators and 
their derivatives at \x\ = L/2. This allows us to express the solution for the n-th transverse 
eigenstate in terms of the field operator ipb at the boundary. Since the results are very 
similar for both contacts we concentrate on the left lead, and the result for the right one 
will be given without derivation. For the contact at x = —L/2 we obtain 

ip{x) = ipbcoskx + jd x ipbsm kx. (11) 

This expression contains both incident and outgoing waves. According to the causality 
principle, the incident wave ipinix) is determined by a state of the lead far away from the 
barrier. Therefore, ipi n (x) must not depend on properties of the barrier. Equating the 



incident part of Eq. (11 ) to the form describing free particles we find 



i An k 2 

where k\ = J2m(e — e n ) and c„ jSi fc is an annihilation operator of an electron in the lead with 
a longitudinal momentum k and spin s in the n-th transverse mode of the lead. 



Eq. (12) relates the field operator at the boundary to the equilibrium states of the n-th 
transverse mode in the lead. We are interested in finding a relation between the boundary 
value of the field operator corresponding to the lowest transverse eigenstate of the conducting 



wire and the incident state of the lead. To find this relation, we project Eq. (12) onto the 
eigenstates of the wire. Since transverse states of the lead are not eigenstates of the wire, 
we obtain an infinite system of linear equations for boundary values of the field operators 
ipj describing different transverse eigenstates j of the wire 

^j-Y, r jj' d ^3' = Zj, ( 13 ) 

r 

where coefficients Tjji depend on the properties of the contact, and we do not calculate 



them as we do not need their explicit form. We must find a solution of Eqs. (13 12) for the 
state j = describing the lowest subband which is responsible for an electronic transport 
in the wire, while the states j > with higher transverse energies do not contribute to the 

10 



transport. It follows from Eq. (13) that the relation we are looking for has a form 



A{e)$ + B(e)dJ = ^= ]T 7 (fc)c n 27r<J(e - e n ), (14) 

V» n=n,fc>0 



where the exact expressions for the coefficients in Eq. (14) depend on the shape of the 
contacts. The boundary condition for the right contact has the same form but with complex- 
conjugate coefficients. 



Coefficients in Eq. (14) are not arbitrary. In particular, they must provide correct anti- 
commutation relations for electronic field operators. It is worthwhile to relate the coefficients 
to such physical parameters of the system as transmission probability t of incident electrons 
for non-interacting system. Therefore, we consider non-interacting electrons, for which we 
can easily solve the equations for the field operators inside the wire. Then we impose a re- 
quirement of fulfillment of anticommutation relations, and calculate the conductance. This 
allows us to reduce the number of undetermined constants. As it is more convenient to 



express boundary conditions in terms of physical values, we multiply Eq. ( 14 ) on the left by 
its Hermitian conjugate, then we transform the obtained equation to the time representation 
assuming that the coefficients are slowly varying functions of energy in the region close to 
the Fermi energy, and neglecting terms corresponding to transitions to upper empty energy 
bands of the wire. Finally, we find boundary condition for the left (right) contact for each 
spin direction 

YP±J + vfIpf = \ £ c+,c n e^'-^, (15) 

n,n' 

where j and p are operators of current and of the smooth part of charge density perturba- 
tions, pF is the 2fcp-component of charge density, which is related to the Friedel oscillations, 
/ is a number of the order unity if the transmission probability is not close to unity, and 



/ ~ v/2(l — t) if 1 — t <C 1. Thus the Friedel oscillations disappear if the contacts are 
adiabatic. 



In order to check the validity of conditions ( 15 ), we considered a wire with non-interacting 



ID electrons which is attached to smoothly widening leads, so that the contacts are nearly 
adiabatic. We also assumed that there might be a potential step of the height Uq <C Ef at 
the interface. Under these assumptions we were able to use the quasiclassical approximation 
in the lead and match the quasiclassical solution outside the ID conductor with the exact 



solution inside the channel. We found that the condition (15) yielded a correct result for 

11 



the conductance G = tGo in agreement with the Landauer formula, where Go — e 2 /h is the 
conductance quantum. 

As we need the boundary conditions for the LL described in the bosonic representation 



we have to present (15) in the bosonised form. We should note that since one assumes a 
linear dispersion of electrons within the LL approach, the theory is valid provided that all 
energies are small in comparison with the Fermi energy. However, the amplitude of the term 
v F fp F which is responsible for the Friedel oscillations is of the order of the Fermi energy if 
/ = 2\/l — t is not small. Therefore, we consider nearly adiabatic contacts with yl — t <C 1, 
and neglect terms of the higher order in /. 

Then we transform in a standard way^ the fermionic operators to charge and spin density 
variables for right (r = +1) and left (r = —1) moving electrons 

Pr= --—7= [Pr,t + Pr,i\ > °V = ~JR Kt ~ PrJ 

and obtain the boundary conditions for bosonic field $ p at the left and right contacts 

v F d x $ p T d t $ p + V2fe F sin(\/2$ p =f k F L) cos V2$ a = P^ R (16) 

v F dJf a T dA + V2fe F cos(V2$ p t k F L) sin y/2$ a = P^ R , (17) 

where P^ ,R = 2nv F N^' R , N^ ,R is the operators of excess number of charge (y = p) and spin 
(v = a) densities in the left (L) and right (R) leads, respectively. The expectation values of 
the operators P^' R and correlation functions of their fluctuating parts 5P^ )R = P^ ,R —(P^' R ) 



can be calculated easily from the right-hand part of Eq. ( 15 ). The average of P^ R for charge 
channel is proportional to the potential JJ^ R applied to the left (right) contact and, similarly, 
the expectation values (N^ ,R ) are equal to the excess spin densities in the leads, hence the 
difference of the expectation values (P R — P%) can be related to "spin bias" V a : 

(Pp' R ) = U^' R /V2, (P^ R ) = UZ> R /V2, 

Correlation functions are identical for both channels and for both contacts, while correla- 
tions between left and right contacts and between charge and spin operators are absent. In 
the frequency representation correlation functions for anti-commutators and commutators 
read 

({6P(u), 6P(u/)}) = Att 2 u coth 4U(^ + u 1 ), ([<SP(u), SP(u/)]) = Att 2 uj5{uj + u'). (18) 

12 



If there is a metallic gate near the quantum wire we must take into account screening by 
the gate. Following the approach of Ref. |2H we find that the screening by the gate results 



in a modification of the factor in the first term of Eq. (16). Then the boundary conditions 
for the case of short-range interaction acquire the form 

^d x % T d t $ p + V2fe F sin(V2<l p T k F L) cos V2$ a = P^ R . (19) 

p 

This modification of the factor before the spatial derivative can be also illustrated by means 

of the simple model in which the factor K p is equal to 1 at the non-interacting lead x = 

—L/2 — and step-like reaches its value in the wire at x — —L/2 + 0. Then we integrate 



the equation of motion (10) from x = —L/2 — to x = —L/2 + and assume that $p is a 



smooth function of x we obtain 



d x %{-L/2 - 0) = ±jd x $ p {-L/2 + 0), 

p 



which explain transition from Eq. (16) to (19). 



In case of a wire adiabatically connected to ideal Fermi-liquid reservoirs at x — ±L/2 the 
boundary conditions reduce to 

fe* T d t ) * v {x= ± L/2)=P^ R . (20) 



Thermodynamically averaged boundary condition Eq. 20 for the charge channel is reduced 
to that derived in Ref. [2U 

It looks natural that in case of gated quantum wire the gate screens externally applied 
electric field and the problem is described in terms of boundary conditions, as it was discussed 
in Ref. EH Of course, inside the wire there is also an electric field induced by non-uniform 
distribution of electrons, but this electric field is taken into account by the interaction 
between electrons. However, it looks less clear whether one can describe the driving voltage 
by boundary conditions when there is no gate (the case of long-range interaction). Therefore, 
in case of long-range Coulomb interaction we considered two approaches. First, we derived 
equations of motion for the phase fields with driving dc voltage taken into account by 
boundary conditions. Second, we inserted the driving dc electric field into the Hamiltonian, 
when the external field appears in the equation of motion for the displacement field & p (x, t). 
But the equation of motion for the displacement field $ p (t) at the impurity site turned 
out to be the same in both cases. For the external ac electric field the result of these two 
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approaches is different, and this is not surprising, because the ac-fieid can be induced by 
incident electromagnetic wave, and this could lead to results different from the case of the 
ac voltage applied to the contacts. 

The fluctuating part of the boundary conditions takes into account that the leads play a 
role of a heat bath. Thus it describes relaxation in quantum wires induced by interaction 
with bulk leads. Such mechanism of relaxation is dominant in practically important case of 
the quantum wire which is shorter than the relaxation length related to coupling of electrons 
to phonons or to other dissipative bosonic bath. 

Let us shortly discuss the role of damping of the plasmonic excitations. The equation 



of motion (10) does not contain damping. If one takes into account relaxation inside the 
wire then a finite damping must appear. Damping due to electrostatic coupling to metallic 
gate was studied in Ref. [31] and its effect on electronic transport passing an impurity was 
considered in Ref. [32j In Ref. [32] the effective action was reduced to a simple form that 
allows to find the Green's function of the system with finite damping. In this case the 
Fourier transformed retarded Green's function acquires the form 

Dg = "' f 22 (21) 

lo z + iujv — q A v l p 

where v is the damping constant. According to Ref. EH coupling to fluctuations in the gate 
is relevant only at K p < 1/2. However, following Ref. [321 one can consider v > at all values 
of K p as a phenomenological parameter bearing in mind that other sources of dissipation 
are possible. Below we concentrate on the case of not very long conducting channel, L <C l u , 
that is on the case which is opposite to that considered in our letter^. Then the damping 
constant v can be neglected. However, even under condition v = some dissipation in 



the wire is still present since the boundary conditions (20) with (18) connect the wire to a 



thermal bath. This can be demonstrated by calculation of the eigenmodes in a pure wire 



without impurity. From Eq. (10) with boundary conditions (20) we find the modes with 



spectrum uj n = q n v p and eigenfrequencies given by relation 

tiv p iv p 1 + K 2 

Similarly, damping of the order of Vp/L is present in case of the long-range interaction. We 
see that damping of the eigenmodes is finite, and the damping constant decreases when the 
length of the wire increases. Further, distribution function of the excitations in the wire is 
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determined by correlation functions ( 18 ) so that in the equilibrium state it is described by 



the Planck distribution. In the non-equilibrium state the excess energy can be dissipated in 
massive leads and if the wire is short enough this restricts heating effects. 

C. Equations of motion of the displacement field at the impurity site 

In this section we consider ideal adiabatic contacts and derive equations of motion for 
the phases $ p and $0- at the impurity. Consider first the case of short-range interaction. 



In order to solve equation of motion (10) for the phase field we use Fourier transformation 



with respect to time and formally solve differential equations for <& v {u,x). Solutions are 
linear combinations of exponents exp (v^j and exp ( — ^p j with different coefficients to the 
left and to the right of impurity. We match the solutions at the impurity site and use the 



boundary conditions (20). In this way we express operators & u (u, x) in terms of their values 
at the impurity site, x = 0, and after inverse Fourier transformation obtain equations of 
motion for the displacement field at the impurity site. The equations read 



d t $ p + —^Z <g> sin V2$ p cos V2$ CT = F <g> P p 

a/2 

W i-» r-~ 

d t <$> a + -= sin V2&V cos V2$ P = P a (t- t L ). 

a/2 



(22) 
(23) 



Here £g> means convolution in time, P a = P^ — P&, and Z(t) and F(t) are defined by means 
of Fourier components 



Z(u) 



K, 



1 — iK p tan ujt L 
' K p — i tan ut L 



F(u) 



K n 



K p cos ut L — i sin Lot L 



LKa 

2vp 



(24) 



These equations resemble equations of motion of two coupled overdamped pendulums with 
some retardation function in the gravity term. Oscillatory dependence of Z(u) and F(u) 
describes multiple reflections of the bosonic excitations of the LL from contacts. This state- 
ment can be illustrated by the expression for Z in time representation 



Zit) = K p 



5{t) + 2Y,r m 5(t- mt L ) 



m=l 



K n 



K n 



(25) 



where r is the reflection coefficient of plasma excitations from the boundary between re- 



gions with different interaction constants. Eq. (25) is symmetric with respect to left and 



right contacts because we consider the impurity located in the middle of the conducting 
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channel. Generalization to the case when impurity position is shifted from the middle is 
straightforward and we do not consider this case. 

Consider now the case of long-range Coulomb interaction between the electrons. If we 
consider interaction potential in the system of finite length rt5J) formally at all values of 
coordinates, then we find that it is symmetric with respect to the left contact and periodic 
with period 2L. This means that we can expand the field operators in Fourier series in the 
interval (-L/2, L/2) 

®n{x) = -° + ^ <l n cos q n (x + -) . (26) 

Z n=0 ^ l ' 

This result in simple and easily soluble equation of motion 

(u 2 -q 2 n v 2 n )K -v 2 F l[dM0) ~ (-l) n dML)] = < ^^W{smV2$ p (t,0)cosV2$ a (t,0)} u . 

(27) 
Then using the boundary conditions we derive equations of motion for the field operators 



at the impurity position x — 0, and find that they are similar to (22 23) but with different 
memory functions F and Z 

iu>R+-2u 2 {R\-R 2 _) _ iuR- v F « (±l) fc 

{) " TT2ujR~ + ' F{U) ~ l + 2zciV R±{U) ~ T A ^ - qW{q 2k ) 



k=—oo 

with q n = j^, v 2 = v F [l + 2 / ~fK Q (\qd\)], where McDonald function describes Fourier transform 
of Coulomb interaction in a ID conductor of transverse dimension d. We will concentrate on 
the case of long enough wires, and we will need functions Z and F at u = 0, when Z(0) = 
F(0) = 1/2, and at u ^> Vp/L, when summation can be approximated by integration. In 
the latter case we find with logarithmic accuracy 

Z(u) « , — , F{u) « — ?^==, Q. = - ! = - (28) 



where g w is a solution of equation to = qu]V(q u ), and 7 is specified by Eq. Q and determines 

to the strength of the long-range Coulomb interaction. 

And, finally, let us discuss equations of motion for the single-mode LL, that is spinless or, 

in other words, spin-polarized LL. This is the simplest case, and derivation of the equations is 

similar to the spinful case. For the spinless LL we concentrate on the short-range interaction 

with constant value of K p . Then the equation of motion for the phase at the impurity site 

reads 

dt$(t)+ dt 1 W i Z(t 1 )sm2$(t-t 1 ) = - dtiF(t-ti)P(ti). (29) 

Jo 2 jo 
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Equations (22 23) and (29) are non-linear equations for operators, therefore, it is not 



easy to solve them in general case. As (29) and (22 23) resemble equation of motion of an 



overdamped pendulum and two coupled pendulums correspondingly, one can expect that 
when the system is driven by a constant external bias the phase can increase and oscillate, 
which means presence of both dc and ac current in the system. The simplest case is when 
the inter-electronic interaction is strong enough (K p <C 1 or 7 is of the order or larger then 
unity). It will be shown that under these conditions fluctuations of the phase field $ p are 
relatively small and can be described by Gaussian approximation. We will solve the problem 
for the case of spinless LL in Gaussian approximation in the next section, and at the end 
of the section we will discuss the validity of Gaussian approximation. The spinful case is 
more complicated as fluctuations in the spin channel at K a rj 1 are not small and are not 
Gaussian, however, for the considered case of spin- independent interaction, K a = 1, they 
can be taken into account strictly by means of refermionization method. 



III. DYNAMIC REGIME OF CONDUCTION IN THE SPINLESS LUTTINGER 
LIQUID 



A. Gaussian approximation 



In this section we will consider the most technically simple case of the single-mode LL 
with short-range interaction between electrons. 

First, we represent the bosonic field operator at the impurity site as a sum of its expec- 
tation value and fluctuating part, $ = $ + 0, $ = ($). Then we perform thermodynamic 



averaging of both sides of Eq. ( 29 ) and obtain equation for expectation value $ of the field 
operator at the impurity site 



d t $(t) + / dhWiZ(t - £i)(sin2$(ti)) 
Jo 



dtiF(t - ti)V{ti 



(30) 



Eq. (30) is not a closed equation for $(£) since it contains an expectation value of sin 2$ (t) 
which depends both on expectation value $ and on fluctuations <p of the displacement 
field. Therefore, in order to calculate the expectation value we need to study fluctuations. 



Subtracting Eq. (30) from Eq. (29) we find an equation of motion for the fluctuating part 
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of the displacement field at the impurity site 

d t <p{t)+ rftiW / i Z(t-ti)[sin2<l>(ti)- (sin2$(ti))l = / dhF(t -ti)SP(ti), (31) 

JO JO 

Equations to be solved somewhat resemble equations for the superconducting phase in 
the Josephson junctions. However, there is an important difference from Josephson junc- 
tions since in the latter case fluctuations of the phase are, typically, small and one can 
develop equations in powers of small fluctuating phase and consider the Gaussian model 
for fluctuations. In our case fluctuations are not small, nevertheless, in order to solve the 
problem we simplify the problem and assume that fluctuations are Gaussian. Strictly speak- 
ing, the fluctuations are not Gaussian, and in general case this is just a model assumption. 
However, we will demonstrate below that this approach can be justified in case of strong 
inter-electronic repulsion and in the limit of high voltages, where the non-Gaussian part of 
fluctuations becomes small. 

Thus we solve the problem by means of the self-consistent harmonic approximation^, in 
which fluctuations are assumed to be Gaussian. In this approximation, we replace 

sin 20^2/i0, h = e~ 2{ ^\ (32) 



and instead of (30) we obtain more simple equation for the expectation value of the field 
operator 

/•oo /*oo 

d t $(t)+ dt 1 W i Z(t-ti)h(ti)am2$(t 1 ) = dtiFti-txWih), (33) 

Jo Jo 

and a linear equation for fluctuations 

~ /"OO A /"OO ^ 

d t 4>(t) + 2Wi dtiZ(t-t 1 )h(t 1 )cos2$(t 1 )4>(ti)= dhF(t -£i)<JP(ti). (34) 

Jo Jo 

Coefficients of this equation depend both on the mean square fluctuations ((f) 2 (t)) and on 



the expectation value $, so it must be solved self-consistently with Eq. (33). 



B. Low applied voltage, V < Vr 



If the applied dc voltage is small enough, V < Vp, Eqs. (33) and (34) have stationary 
solutions for phase $ and for mean square fluctuations (<f) 2 )- In this subsection we consider 



this stationary solution and find Vy. In the stationary case Eq. (33) reads 



Wihsm2<5> = V, (35) 
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and Fourier transformed Eq. (34) reduces to the simple form 

- iuj)(tu) + 2WihZ(co) cos2$0(w) = F(u)5P(u). (36) 

This equation can be solved easily. Then taking into account correlation functions given by 



Eq. (18) we can calculate mean square fluctuations 



Using Eqs. (36) and (24) the expression under the integral sign can be presented in a 



convenient form 

(fc\ = I Sl [°° cjcothj^o; 

K9) 2 Loo(u* + W?) [{l + KD + {l -K$)sm(ut L - a u )]' { ] 

where 

W 2 -u 2 
a^ = arctan — p— — — , W c = 2WiK p h cos 2$. 

2uW c 

The result of integration depends on relation between Vr and temperature T. First, we 
consider the limit of zero temperature. In pure LL this integral would diverge logarithmically 
both at high and low frequencies. The divergence at the upper limit in the TL formalism 
must be cut off at frequency A of the order of the bandwidth or the Fermi energy. The 
infrared divergence at low frequencies is a distinctive feature of ID systems. We find that in 
the presence of impurity the infrared divergence is cut off at frequency W c ex W». In addition 
to the logarithmically divergent part, the denominator contains the oscillating factor induced 
by reflections of fluctuations from contacts. If the length of the quantum wire is large enough, 
L ^> v/Wih, the main contribution to the integral is determined by frequencies uti ^> 1 
and oscillations contribute little to the integral. Then the oscillating factor can be replaced 
with its average value K p /2, and integration yields 



Since W c depends on (0 2 ), Eq. (38) determines the self-consistency condition for (0 2 ). Sub- 



stituting h from Eq. (32), we find 



K A 

1 2(1 -K p ) 2W l K p cos2& v ; 



Now using this equation we can calculate maximum value of the left hand side of Eq. (36) 



which determines the value of the threshold voltage below which the static solutions for mean 
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phase $ exist. We find that H^/tsin2$ reaches its maximum value when cos 2$ = jK p , 
and expression for Vt has a form 



Vt = 2W t I — ^ j y/l - K p . (40) 

We see that the threshold voltage at low temperatures is determined by the impurity po- 
tential renormalized by quantum fluctuations. In case of interelectronic repulsion, K p < 1, 
the mean square fluctuations ((f) 2 ) and, hence, Vp are finite, while in non-interacting system, 
when K p = 1, fluctuations become infinite and Vp is destroyed by quantum fluctuations. 
Thus we find that the solution for $ is stationary at V < Vt, that is current cannot pass 
an impurity. This result is a consequence of our approximation in which only Gaussian 
fluctuations were taken into account. If we took into account fluctuations of solitonic type 
for which the phase increases by tc due to tunneling, we would obtain a small tunneling 
current at V < Vt- Fluctuations of this type yield the well-known power-law I-V curves 1 . 

Now we consider the case of finite temperatures. The self-consistency equation has so- 
lutions which correspond to a finite value of fluctuations only if T < T c ~ Vp,o = Vp(T = 
0, L — oo), so there is a characteristic temperature above which Vt is destroyed by thermal 
fluctuations and the impurity does not suppress electronic transport. 



If the quantum wire is short enough, L ~ v p /Vr,o, we must not average Eq. (37) over 



oscillations at utp ~ 1. At these frequencies {(p 2 (u)) in Eq. (37) is proportional to u;" 1 as 
before but with a different factor. As a consequence Vp is suppressed in short wires, and 
impurities do not destroy the linear conduction when L < L c ~ v/Vpfi- This happens due 
to increase of fluctuations at the impurity site because fluctuations are reflected back from 
the contacts, while the distance to the contacts becomes smaller than the correlation length 
of the fluctuations. 

Thus we have obtained that conductance is greatly suppressed below the threshold field 
Vt which is determined by impurity potential renormalized by fluctuations. As fluctuations 
depend on temperature and length of the wire, Vp decreases as temperature increases or the 
wire becomes shorter. Results of numerical calculations of Vt as a function of temperature 
and wire length for K p = 0.3 are presented in Fig. fl] Results for other values of K p are 
qualitatively similar. One can see that there are critical temperature and the inverse wire 
length above which there are no solutions with non-zero Vt- 
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FIG. 1. Results of numerical calculations of the normalized threshold field as function of temper- 
ature and wire length 

In Fig. [2]we present a "phase diagram" which shows characteristic temperatures and wire 
lengths limiting the region Vr > where conductivity is suppressed by impurity. 

C. Dynamical regime near the threshold voltage 



Now we consider the behavior of the system at voltages slightly above the threshold 



value at low temperatures. As equation of motion (33) resembles equation of motion of 



an overdamped pendulum, near the threshold the phase must increase very non-uniformly, 
since under such conditions the pendulum would slowly move up and rapidly fall down. So 
we expect that when 5V = V — Vr <C Vr, the phase increases very slowly if sin 2$ is close 
to unity, and rapidly passes the values of $ for which V — Vr sin 2$ is large in comparison 
with SV. Then as the phase increases slowly, fluctuations can be described by means of 
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FIG. 2. Region of temperatures and wire lengths in which the impurity strongly affects the elec- 
tronic transport 



quasi-static adiabatic approach and we can use Eq. (39). The latter is applicable provided 



the values of $ satisfy condition cos 2$ > 0. Note that the maximum value of Wj/isin2$ is 
determined by cos2$ m = jKpi that is cos 2$ is positive in the vicinity of the maximum. 
Then in the region around the maximum where the phase spends the most part of the period 
we can, indeed, use the adiabatic approach. In this region we can develop Vt sin 2$ near its 
maximum in series with respect to 5$ = 2($om — 3>) and obtain equation of motion for the 
phase in this region 

d t 5$ + V T ( 1 - — ] = V. (41) 



\-K n 



Its solution around the maximum has the form 



i(V-V T )(l-K p 



V T 



tan 



\ 



(V - V T )V T 
l-K n 



t. 



(42) 



1-K„ 



(V-Vt)V t 



and then 5$ increases 



We see that 5$ is small during large time to ~ tt 
rapidly during short time, so that in the first approximation in J{V — Vt)/Vt the period 
of oscillations is determined by time t . Hence, with the same accuracy we can calculate dc 
current 

- 2G„,\ <- V r V j} VT . (43) 



\ 



1 - A'„ 
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We can also conclude that time dependence of the current consists of narrow peaks of height 
~ GqVt and width ~ 1/Vr separated by long time to during which the current is small. 

D. I-V curves and noise spectrum at high voltages 



As it was noted already, it is difficult to obtain I-V curves at low voltages accurately 
because of time dependence of the mean square value of fluctuations. The problem is sim- 
plified at high voltages, V ^> Vp, when the mean square value ((f) 2 ) becomes nearly constant 



with small oscillating component. In this case Eqs. (33 34) can be solved perturbatively 



assuming that the oscillating parts of both mean square fluctuations ((f) 2 ) and of the mean 
phase $ are small. 

In this subsection we consider the limit of relatively long conducting channel, Vpti 3> 1, 
but not too long, so that the wire is short in comparison with the damping length related 
to relaxation due to coupling to phonons etc, L <^l u . In this case we have to use the exact 



form of Z(t) in equation for the expectation value Eqs. (33) but can keep only the first 



delta-function in kernel Z(t) in equation for fluctuations Eq. (34). In time representation 



this means that we take into account current pulses reflected from the contacts but we ignore 
correlations between fluctuations shifted by time nip necessary for the excitation to return 



to the impurity after multiple reflections from the contacts. Then Eq. (34) acquires simple 
form and can be solved easily 



d t <j> + W(t)<j>(t)=f{t), W(t) = 2K p Wih(t)caa2$(t), f(t)= / dt^it -t^dPitJ. 

Jo 

(44) 



Solution of Eq. (44) has the form 



i = f dtl f {tl) e-^ VI "' d \ 



(45) 



Using Eq. (45) we can calculate mean square fluctuations (0 2 ). As we consider the long 
channel we average, again, over oscillatory factor in F(t) and find 

tO - f* W{t 2 )dt 2 -[* W(t 2 )dt2-uj(t x -t3) /,„x 

— e J4 i Jt 3 . (4oJ 



l ) 



— p / dti j dt 3 / duu coth — — e 

4 J—oo J— oo J 2,1 



To solve this equation we need to calculate, first, W(t) which is determined by fluctuations. 



In order to do this we solve Eqs. (33) and (46) for fluctuations seeking for ((f) 2 ) in the form 
((f) 2 ) = ((f) 2 ) t + ccoswo^ + ssincuo^, where u = 2nl, and (■ ■ -) t denotes averaging in time. 
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We assume also that c, s < 1. Substituting this form into Eq. (46) and keeping only leading 
terms we obtain in the limit of low temperatures 



2 > 






4 



\uj\du) 



oo00 2 + b 2 



1- 



2W [6(wo + w 2 ) cosu t + w (w 2 — u 2 ) sinu t] 



Wn 



oo 2 ) 2 + Ab 2 u 2 



(47) 



where W Q = 2WiK p e~ 2 ^\ b = (W(t)) t 
obtain 



|c|Wq. Performing integration in Eq. (47) we 



2 ) 



K n 



A nW 2W u q . 
in cos Wo^ In ~r sin w £ 

W CJq o 



(48) 



Thus we have found that the main logarithmic contribution to (0 2 )t is determined by 



relation similar to Eq. (38) valid in case of small voltages, but with different infrared cut-off 



frequency b which is much smaller than W c in Eq. (38). From the self-consistency condition 
we find 



7iK p W 
2u 



K P W ^ 2u 2 



U) 



nW 2 ' 



W Q = V T 



K n 



1-K„ 



-JKf/ T 
2/1 ^K~ P V 



< V T 



(49) 



Here we have expressed Wq from Eq. ( 40 ) in terms of Vr at zero temperature in the limit of 
the long wire. 

We see that at high voltages the solution with finite amplitude of the oscillations exists 
only at K p < 1/2, i.e., when inter-electronic interaction is strong enough. The result differs 
from that for the regime of small voltages when fluctuations do not destroy the dynamic 
regime at any repulsion strength, K p < 1. The result also differs from the case of dynamic 
regime of conduction induced by Friedel oscillations at sharp contacts to quantum wire®!, 
where the critical value at high voltages is K p = 1/3. 

So far we considered in this subsection long channels at the zero temperature limit. As 
fluctuations at high voltages become stronger, they decrease the "critical" temperature and 
increase the "critical" length in comparison with the case of small voltages considered in 



subsection |IHB| Therefore, in the limit of high voltage and low temperatures we obtain 
T c ~ W and L c ~ vV/W 2 . 



Now, when we have found mean square fluctuations (48), we can solve Eq. (33) pertur- 



batively in the limit of high voltages, V 3> Vr, and calculate current using Eq. (pO). The 
total current passing the impurity consists of dc part, I = VG$ + I n i, where I n i is non- 
linear correction to Ohm's law, and of ac part, 7 ac sino;ot, which oscillates with frequency 
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wo = 27r//e ~ eV/Ti (in dimensional units) 



GnW< 



o^o-t 



{l + K 2 ) + (l-K 2 p )cosu t L 
\(l + Kj)-(l-Kl) C osu Q t L 



2G n W 2 



i- ni 



O^o 



V 



\ 2V 2 

111 ^wl + 



(50) 
(51) 



The oscillating factors in these expressions are due to reflections of current pulses generated 
at the impurity from the contacts. Non-linear correction to Ohm's law is shown in Fig. [3} 




6.0 ^ 



FIG. 3. Dependence of non-linear part of the current, I n i, on voltage and interaction parameter 



The ac current depends on the coordinate along the wire, and its value at the contacts 



to the wire is different from that at the impurity site given by Eq. (50). The amplitude of 
the ac current at the contacts reads 



/. 



V2G W 



ac— contact 



l + K 2 -(l-K 2 )cosu t L 
Now we calculate the noise spectrum around the oscillation frequency 



(6I(u)6I(u/)) = dt dt'(SI(t)8I(t'))e 



jj\\ „jajt 



dcu' 

— < 
2vr 



iuj'f 



Then we express SI(t) = dt^/ir from Eq. (31) and find that 



{Si(t)SI(t')} = (W / i /7r) 2 [(sin2<l(t)sin2<l(t / )) - (sin 2$)' 



(52) 



(53) 
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We calculate correlation function in the right-hand side of Eq. ( 53 ) assuming that fluctua- 
tions are Gaussian 

(sin 20(£) sin 20(f)) = \ cos2($(t) - $(£')) e 2 ^«-^')] 2 > . (54) 



Correlation function of the phases can be found from the first term in Eq. (47) that deter- 
mines the leading contribution 

(0(£)0(t')) = -^ [e' bT Ei(bT) + e bT Ei(-bT)] , r = t-t'. (55) 



Using Eqs. (52 55) and calculating the integrals we find two maxima of the noise spectrum 



around frequencies uj = ±u 

isft \sft r\\ ^2-tMi-Kp) itT{l-2K fJ )am.T{K p 

(6I(u)6I(u )) « G 2 V T ( p> — p V^T^ + u) - (56) 

2(1 — K p y Kp K p p \\uj\—ujq\ v 2K p 

Note that the maxima are present under the same condition K p < 1/2 for which the solution 



with finite amplitude of the oscillations at high voltages was found. According to Eq. (56) 



the integral noise power is of the order of ~ G^Wf which is much larger than the ac signal 
power ~ I% c at frequency u . 

In case of long-range Coulomb interaction the correlation function can be found similarly 



to derivation of Eq. 55 , it has the form 



- ~ 1 A Jin 3 A 

(0(£)0(£')) = —In In V 



4 7 Wf\t-t'\ 



Then instead of Eq. 56 we find at u ^> Vp/L when Z(u) is described by Eq. (28) 



i 



, r?/ , rf/ A . Wf /, |u;|||u;| -a; |\ 4 " 7 ., /N 

The noise spectum decreases from its maximum value approximately as 1/u. This is some- 



what similar to the case of very strong short-range ineraction [K p — > 0) in Eq. 56 



E. Validity of Gaussian approximation 

Now we discuss conditions under which the Gaussian model that we have used to describe 
fluctuations can be justified quantitatively. Note that fluctuations of the displacement field 
in pure ID system are Gaussian because the TL Hamiltonian is quadratic, and mean square 
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fluctuations are infinite, (0 2 ) = oo. Impurity makes fluctuations at the impurity site finite 



(confer Eq. 38), but fluctuations become non-Gaussian because of the cosine impurity term 



in the Hamiltonian. As the current passes the impurity, the impurity term oscillates, and 
frequency of the oscillations increases with voltage increasing. This results in a decrease of 
the time-averaged impurity potential making the impact of the impurity effectively smaller. 
Therefore one should expect that relative contribution of the non-Gaussian part to fluctu- 
ations must decrease in comparison with the Gaussian part. Then at voltages 7 > Vy we 
can try to calculate non-Gaussian contribution to fluctuations perturbatively. 

We select two contributions of the fluctuating part of the phase, <j> = (f>a + <f>i, where the 



first term is the Gaussian contribution which satisfies simplified Eq. (34), while (p satisfies 



full Eq. (31). Considering non-Gaussian part <j>\ as a small correction we linearize Eq. (31) 
and obtain equation for <pi. We consider, again, zero temperature and long conducting 
channel, Vt^l ^> 1, when Z(t) ~ K p 5(t). Then equation for non-Gaussian corrections reads 

-wfa + W(t)fc(t) = WiK p x 

{cos 2$(£) [2h(t)j> G (t) - sin 20 G (£)] + sin 2$(£) [(cos 20 G (£)) - cos 20 G (£)] } . (57) 

Next, we derive equation of motion for the third cumulant C3. In the first approximation, 
Cz(t) = (0i(£)0g(O) 2 ). We multiply Eq. (57) by G (£ = 0), and make thermodynamic 



averaging using for the right-hand side the rules valid for Gaussian fluctuations. Then we 
obtain 

d t C 3 (t) + W(t)C 3 (t) = 41V^ p /i(t)sin2$(t)({0 G (t)0 G (O)}) 2 . 

Solution of this equation has the form 



Cs(t) 



dt 



- f* W{t 2 )dt 2 



4W^/i(t 1 )sin2$(t 1 )({(/. G (t 1 )(/. G (0)}) i 



Calculating the integral, using Eq. (55), and keeping the leading terms we find 

AV 2 



C 3 (0) = a 3 K p 



1-KAn 



p ^nWZ 



a 3 = — / dxe~ 
4n Jo 



e- x Ei(x) + e x EU-x) 



0.3490658. 

(58) 



Similarly, from Eq. (57) we can calculate the fourth cumulant C 4 = (0i0 G (O) 3 ). The 



result for the cumulant at coinciding times is 

1 r° 



C 4 (0) = a 4 K A a 4 



16 



e~ x Ei{x) + e x Ei(-:c) 



e x dx w 2.7224372. 



(59) 
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Now we can compare non-Gaussian contributions (58, 59) with related Gaussian con 



tributions. The third cumulant contains two terms. The first one does not increase with 
voltage, while the second term increases with voltage but more slowly than (0 2 ) given by 



Eq. (48) and contains additional small factor K p . Non-Gaussian contribution (59) con 



tains both small factor K p and does not increase with voltage. Thus we see that non- 



Gaussian contributions (58, 59) are relatively small. For example, the non-Gaussian con- 
tribution to (0 4 ) given by Eq. (59) is small in comparison with Gaussian contribution, 

3(4>r~K^ 2 (i)G 



V \ \-2K p 

V T 



which increases with voltage. 



IV. DYNAMIC REGIME OF CONDUCTION IN THE SPINFUL LUTTINGER 
LIQUID 

A. Refermionization in the spin channel 



In Sec. Ill we treated fluctuations in the Gaussian approximation which can be justified 
only in the limits of strong interaction at high voltages. In the spin channel, fluctuations 
at the impurity site are always strong and non-Gaussian. Therefore, though the Gaussian 
approximation can be used as a model one, it is not justified. However, if interaction is 
spin-rotation invariant (K a = 1) and the impurity is situated in the middle of the wire 
we can make use of the symmetry of the system and solve the problem strictly by means 
of refermionization method. This method consists in introducing new fermionic variables 
for spin channel. Equations of motion for these variables are linear, and, hence, soluble. 
Refermionization was used successfully to treat charge fluctuations in the spinless case for 
the specific value of interaction parameter K p = 1/2 (Ref. [2H and [35]) and to describe spin 
fluctuations in the spinful case for K a = 1 (Ref. ESI and EH|) ■ Following the approach of 
Ref J2H we introduce new phase fields by relation that differs from the similar one in Ref. 
by a constant factor 



Mx 



V2 



®Jx) + eJx] 



± 



i 



V2 



$, 



-x) - QJ-x) 



(60) 



New fields are completely decoupled and the impurity term couples to the field 0+ only. 
Then following Matveev^ and Egger and GrabertP^l we introduce new fermion variables 
v 2™ e "^ + = 9^1 9 = ^ + $ i where g/y/2 is an auxiliary Majorana fermion operator. Then 
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we derive equations of motion for Heisenberg operators if) and find that they depend on 
x — vpt. Equations of motion for operators ipi t2 {t) = ip(x = =F0,£) at the impurity site and 
for g have the form 

vf(4>2 ~ ipi) = igf, vpi^l - tpl) = -igf (61) 

d t g = i{f(4> 1 +4> 2 )-f(ti + 4)}, (62) 

where f(t) = y/2-naW cos V2& p . 

Density perturbations of new fermions are related in a standard way to the gradient of 
the displacement field 

it$* - {4>2$x)o = 7r d ^+( x )- ( 63 ) 

We will consider the limit of strong interaction between electrons when fluctuations in 
the charge channel are small. We represent the field operator at the impurity site again 
as a sum of its expectation value and fluctuating part, $ p = $ p + <p p , $ p = ($ p ) and take 
into account fluctuations <j) p only up to the first approximation. Then commutators of / at 



different times are small and we can ignore time-ordering and solve equation (62) for g(t) as 
if f(t) is a c- number 

g = 2i rdti[/(ti)A(ti)-/(ti)^I(ti)]exp[-— f f{t 2 fdt 2 ]. (64) 



v F Jt 



Now by means of Eq. (61) we can express operators ip(x = 0) = {ip\ + il>2}/2 ■ Then using 



e^+ = V2iragi>= -V^g^ + ?p 2 ), e~^ + = -\[2xa0 = — y/2xag$\ + $) (65) 



we can find cos y/2<& a which we insert into the equation of motion for the charge phase (22) 



cos V2$ ff = J 7 ^ g(4> - 4>l) (66) 



2 

where we made use of the anticommutator relation {g(t), V4(i)} = — • 
We obtain the following expression for cosv2$ CT : 



cosy/^&vit) = 2maW f dt x cos V2® p (h) (67) 

J — oo 

x {[Mt,) - $(*i)]&(*) + $(i)$i(*i) - $(*i)]} exp [- — f f(t 2 ) 2 dt 2 }. 



We insert ( 67 ) into the equation of motion for the charge phase ( 22 ) . We consider the limit 



of small fluctuations. In this case averaging over charge and spin variables can be performed 
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separately since the fluctuations in spin and charge sectors are independent. Expectation 



values of fermionic densities in (67) can be associated with distribution function of new 
fermions by relations 



',t 



imiiMh)) 



ds 
2tcvi 



-n(e,t)e l ^- t2 \ (^fa), #(*„)> 



ds 

2-KVi 



[l-n(s,t)}e~ t£{tl - t2) , 
(68) 



where t = (ti+t2)/2, and other pairings are equal to zero. This is true because operators with 
subscript 1 are related to the incident spin excitations. These excitations are not affected 
by the impurity because the coefficient of reflection from the contact r = \r^f- is equal to 
zero for K a = 1. Note that this case differs from the problem for charge channel considered 
in Ref. |2H Charge excitations incident on the impurity contain fraction of excitations which 
reflected from the contact and, therefore, interacted with the impurity. 

Now we need to determine distribution function n(s,t). To do this we, first, subtract 

-L/2 for spin sector and obtain 



boundary conditions (20) at x — —L/2 and x 

v F d x 4 



L 
-2'* 



P„ 



(69) 



Then we express the derivative d x <p + in terms of fermion density from (63) and take the 



expectation value. In this way we find the condition for the fermion density expressed in 
terms of the distribution function 



ds 

2^ 



n{e,t)-n F {e)] =V tT (t). 



(70) 



Next we multiply equations (69) taken at different times t\ and ti and calculate the expecta- 



tion value. Reducing products of four fermions to sum of products of pairs in a standard way 



and using correlation functions (18) we arrive at the kinetic equation for the distribution 
function 

J n{s - w, t) [1 - n(s, t)]ds = |(l + coth ^) . (71) 



Solution of equations (70 71 ) is 



n(s,t) 



e = Vgjt} 

1 + e t 



(72) 



Thus distribution function of new fermions has a form of the equilibrium function with the 
chemical potential equal to spin bias. Note that the distribution function has such a form 
because at K„ = 1 there are no reflections of excitations from the contacts. In case of 
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spinful electrons with K p = 1/2 we would obtain kinetic equation different from (71) which 



does not have solution in the form of the equilibrium distribution because particles, incident 
on the impurity, contain a fraction that passed the impurity and then reflected from the 
contact. Reflection coefficient in this case is equal to r = ~ K P 
the distribution function of fermions assumed in Ref. ED needs some justification. 



~. Therefore, the form of 



Using equilibrium distribution function (72) in (68) we insert (67) into (22). Then we 



perform integration over energies and find closed equation for the charge phase 

i W _ . /-- ,,. f°° „ TcOS\/2& p (t -ti) -2w f' cos 2 V2$ p (t 2 )dt 2 



h 



x cos V a ( t - ^ ) ti 



F®P„ 



(73) 

where w = 27raW /2 /-UF is the characteristic potential which results from the impurity po- 
tential renormalized by spin fluctuations. This expression is strict in the limit of strong 
interaction, and now we will discuss the conditions of validity of our approach that assumes 
smallness of the fluctuations at the impurity site. 



To estimate fluctuations we will simplify equation (73) taking into account logarithmic 



divergence of the integral at t± = 0. Then with the logarithmic accuracy we perform inte- 
gration neglecting t\ dependence of the regular part of the integrand and using the standard 
ultraviolet cut-off of the integration at t\ ~ 1/A. We obtain 



d t ^ p + V Z <8> sin 2^$ p = F®P P , V 



w . A 
z= In — . 



(74) 



This equation is similar to Eq.(29) for a single- mode LL and can be made identical to (29) 



by changing notations. Therefore, for the case of short-range interaction we can use the 



results of Sec. III. Then using Eq. (38) we find that in the limit of low voltages fluctuations 



are small provided K p hi - C 1, while from Eqs. (48 49) we find that in the limit of large 



voltages fluctuations are small under condition K p \n ^- <^i 1. 

Let us estimate now the strength of fluctuations in case of long-range interactions. We 
solve Eq. 



74 



in linear approximation in fluctuating part of the displacement field 



$ p - <$> p , $ p = ($ p ) and find 



F ^ 6P " C = 2V2V (Z®cos2V2$ p ) t , 

—iu + C 



where ()t means time-averaging. To calculate constant C we must solve Eq. 74 for expecta- 



tion value <3> p . Here we will assume that temperature T is low enough, T ^C Vq, and limit 
our estimation by the cases of small and large voltages. 
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According to study of the dynamics in Sec. Ill, we can conclude that at small voltages 
below the threshold $ p is constant. Then we obtain with logarithmic accuracy again 



<**;> 



-Sin ^ 

4 7 I dV / 



In the limit of large voltages the phase increases linearly 2v / 2$ 
2ttI ?a 2V. Estimating F(oj) and Z(uS) we find 

Vt to 



u t with cjq = 2nf 



C 



(L 



^v 72 ^ 1 ^- W 2aln ^ 

Then using Eq. [18] we obtain with logarithmic accuracy 



-Ilnln'^ ^ 
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From these estimates we conclude that in case of long-range Coulomb interaction fluctuations 
of the displacement field at the impurity site in the charge channel are not large even at 
moderate values of parameter 7 of the order of unity or larger, which correspond to typical 
values of Fermi velocity (confer Eq. Efc, and under this condition equations of motion can 
be treated in quasiclassical limit. Though fluctuations increase with voltage increasing the 
dependence is very weak since ujq is inside the double logarithm. 



B. I-V curves and current oscillations 



To calculate current in the system we must solve Eq. 73 in which we neglect fluctuations 
and substitute $ p by its expectation value <& p . But it is not simple to solve it analytically. 
Therefore, we restrict our study by limiting cases. 

The simplest case is the regime of current bias. We consider this regime and assume that 

V a is time- independent. 

UJq w '' 00 



V(t) = 1 — 1= / at Y(t) sm / dt\ cos cos V a t\ 

2 \/2 Jo 2 Jo 2 



T { 2w uioti 

x , exp <^ -wti sm -—- cos 

smh irlti l loq 2 

where Y(u) = Z(u)/F(u) . 



Uoh 



+ LU (t-i) 



(75) 



If we perform time averaging of (75) we find the static I-V curves. 

Wo 



V dc = - 

4 Jo sinh-TrTtx 



(76) 



smwoti h \ — sm — - + sm — - J — sm — — 



-wti 
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It is worth noting that the result is the same for both short-range and long-range interac- 
tion, which is not very strange since we consider here the limit of strong interaction, and 
fluctuations in the charge channel are neglected. 

I-V-curves can be presented in simple analytical form in limiting cases. Consider T <Cw. 
In the limit of small current, ujq <C w, the second term dominates over the first linear term 
and current is proportional to J—- This corresponds to the result I ex yV expected for 
the case K a = 1 and K p — > 0. 

The results for the opposite limit of large currents, Uq 3> w, are similar in cases of both 
voltage and current bias. At low temperatures, T <C w, we obtain 



TA U ° , W 



1 4uqW 

— arctan ■ 



IX 



ojI - V} 



(77) 



7 o 

The asymptotic I-V curve is parallel to the Ohm's law corresponding to conductance quan- 
tum 2G = e 2 /{nh) with the excess voltage V exc = |?. The excess voltage does not depend 
on temperature, at temperatures T > w its asymptotic behavior is 

v exc = -^ tanh ^- 

The general view of the I-V curves is presented in Figs. |4} 



Time dependence of voltage ( 75 ) can be characterized by amplitudes of harmonics n > 
which for V„ = are 



Vn 






2 4 Jo sinli7rT£i 

1 r 1 i 



Sinw tl I\n+\{z) + I\n-l{z) + - COS C^l I\ n +l{z) - I\ n -l{z) \ 



-wti— vn "„ x 
e ) 



with 



2w . ojoti 
z = sin 



u 2 

Here the difference between short-range and long-range interactions is present since Y(w) = 
Z(uj)/F(oj) in these cases is different. 

At small currents, uq <C w, amplitude of harmonics decays slowly, approximately as 
1/y/n. In the limit of large currents, Uq ^> w, harmonics decay as power law, and with 
logarithmic accuracy we obtain for n > 

87T \2u J UJq 
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FIG. 4. I-V curves at different temperatures and spin bias. Voltage is measured in units of the 
characteristic potential w, and current in units of w/Gq. Dotted line: Ohm's law / = 2GqV. Solid 
lines: V a = and different temperatures T = 0, 0.2u;, l.Ow; from bottom to top. Dashed lines: 
T = and different values of spin bias V a = O.lw, 0.3w, 0.5w from bottom to top. The initial part 
of the I-V curves is shown in the inset 



Consider now the case of the voltage bias when the system is driven by external voltage 
V + V\ cos ut, and assume the limit of large voltages, V ~ ujq ^> w. In this case we can 



consider the second term in the left-hand-side of Eq. 73 as a small perturbation. If the ac 



voltage is absent we obtain the result coinciding with (77). Small ac voltage modifies I-V 



curves, and the most impressive part of this modification is the resonant steps analogous 
to the Shapiro steps in the Josephson junctions 37 . In contrast to Josephson junctions these 
steps are not at constant voltage, but at constant current I = ef like in the regime of 
Coulomb blockade^ and in the regime of sliding CDW in linear-chain conductors^. At this 
current the frequency of the ac voltage is equal to the frequency of current oscillations in 
the wire. The width of the step at V ^> w and Pi<V can be calculated straightforwardly 
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using perturbative approach. We find with logarithmic accuracy 

VxW lT1 , .. A 

V s tep = -frfWo ln — • 
TTV CUo 

Non-zero spin bias induces a spin current which contains both dc and ac parts. The spin 
current can be calculated according to relation 

71 



where $ CT can be found from equation of motion (23) with sinv^^ which is calculated 



using equations (65) and (68). Then we obtain for spin current 



Vff W fc*. t.\ f°° i. fc*. i, . sTsOlVffti -w f* cos 2 V2® p (t 2 )dt 2 ,„ n . 

Z7T Z7T Jo sinn 7ri ti 

In the limit of large voltages, u> ^> w the spin current can be presented in a simple analytic 

form 

K/, w . \ 
I a = — - H sin wot • 

27T V 7TCJo / 



V. NON-IDEAL CONTACT 

A. Refermionization in the spin channel 

As non-ideal contacts induce Friedel oscillations in the quantum wire, one can expect that 
such contacts must induce the effects in transport which are similar to those in the system 
with impurity studied in previous sections. This statement is supported by the results of our 
letter ^ where we have studied the spinless LL with two identical non-adiabatic contacts. 
However the problem of transport through non-adiabatic contact to a quantum wire with 
a spinful interacting electron gas was not solved. In this section we consider electronic 
transport through a clean quantum wire described as a spinful LL with one ideal adiabatic 
and the second non-ideal contact. The main difficulty in solving this problem is, again, large 
fluctuations of the displacement field $ CT . And, again, we solve this problem by means of 
refermionization in the spin channel. 

To study the role of the non-ideal contacts we act similarly to the previous sections, we 
solve equation of motion for the displacement fields with boundary conditions. 
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Consider boundary conditions for the spin channel with K ff = 1 with ideal adiabatic 
contact at x = L and non-adiabatic contact at x — 0. Then boundary conditions read 

(v F d x - dt)K(x = 0) = P CT L - V2fe F sin V2$ a cos V2$ p (80) 

(v F d x + d t )<5> a (x = L)=P R . 

As field & a (x,t) satisfies the equation of motion 

(v 2 F d 2 x -d*)K(t,x) = 0. (81) 

We can find solution for $ CT (x, t) in terms of its values at the contacts, and using then 



boundary conditions (80 ) we can obtain equation of motion for the displacement field at the 
non-ideal contact 

dA + s/2fe F sin V2$ a cos y/2$ p = ^(t) - F?(t - t L )]. (82) 

This equation resembles equation of motion for the displacement field at the impurity site. 
In order to map the problem of non-ideal contact to the impurity problem we consider the 
model LL with ideal contacts at x — ±L and an impurity characterized by back-scattering 
matrix element W at x = 0. Then equation of motion for the displacement field and 
boundary conditions for such an impurity read 

(v 2 F d 2 x - <9 t 2 ) $ CT (t,x) = V2v F WsinV2$ a cos V2$ p 6(x), (83) 

(v F d x Td t )K(x = TL) = Q L > R . 

Here we denote external sources of fluctuations as Q, and later we will relate them to the 
source terms P. 

Next, similar to the case of contacts we derive equation of motion for the phase at the 
impurity site x = 

d t $ a + -j= sin V2$ a cos V2$ p = 2 \Q L (t - t L ) - Q R (t - t L )\. (84) 



Now we compare equations (82) and (84) and find that equations of motion become 
identical if we choose 

W = 2fe F) Q L (t) = P L (t + t L ), Q R (t)=P?(t) 

Thus using such substitutions we can use the results for spin channel obtained in Sec IV A 



to study the problem of electron transport in quantum wire with one non-ideal and the 
second ideal contacts. 
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B. Solution for the charge channel 



Now let us consider the charge channel. Following the method used in Sec. II C we find 
solution of equation of motion for & p (x,u) satisfying the boundary conditions. In this way 
we obtain expression for the displacement field which depends on the values of both $ p 
and i&a(x,u)) at the boundary with non-adiabatic contact, as both these fields are present 
in the non-linear term of the boundary condition ((171). Then using this solutions at x = 



we find non-linear equations of motion for <£> p (x = 0) which are similar to (22), but with 
different memory function Z and different right-hand containing the source terms Pl,r in a 
non-symmetric way 

d&p + fe F Z <g> sin V2$ a cos y/2$ p = F L <g> P p L - F R ® Pf. (85) 

For the short-range interaction Fourier components of the memory functions read 

Z{u) = F L {u) = K p — — g ^r^———- , F r (uj^ - 



2K p -i{l + K 2 p )t&Yi2ujtL 2K p cos2tut L -i(l + K p 2 )sm2ut L 

In the limit of strong interaction between the electrons (K p ^C 1) these equations will give 
the results similar to the case of impurity. 



In case of long-range Coulomb interaction between the electrons we act as in Sec. II C 

and find 

7( s F ( x^ iuQ + - u 2 (Ql - Q 2 _) iuQ_ 

{U> L{ ' \ + 2iujQ + -uj\Ql-Qiy R{UJ> l+2iuQ + -u 2 (Ql-Q 2 _Y 

O (co) - VF V {±1)k 

In limiting cases u = and u ^> vp/L which is a proper limit for long enough wires, we 
find relations for memory functions similar to the case of the impurity 

r,, s ,-, / n 1 ^ / s ^QujL to 

Z(u) = F L (u) « . — , F R (u) « r== , q. 



where g^ is solution of equation a; = q u v(q w ). 

Thus we find that the problem of electron transport in quantum wire with one non-ideal 
and the second ideal contacts is mapped to the impurity problem. Then all the results 
obtained in Sec. HV] can be used for the case of non-ideal contacts if we substitute the 
impurity potential W for fep which characterizes the amplitude of the Friedel oscillations 
induced by the non-ideal contact. 
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VI. CONCLUSIONS 

Using the approach based on the bosonised Tomonaga-Luttinger Hamiltonian we have 
studied electronic transport in ID conductors with a single isolated impurity or with non- 
ideal contacts to leads of higher dimension, and predicted a new dynamical regime of con- 
duction in which applied voltage in addition to dc-current induces ac oscillations with the 
wash-board frequency / = I/e. This regime is different from the well-known regime of 
macroscopic tunneling leading to the power-law I-V curves and occurs in case of repulsive 
interaction between the electrons and is related to sliding of the Friedel oscillations induced 
by the defect. 

We have considered both short-range interaction and long-range Coulomb interaction. 
The former appears in gated quantum wires where the long-range part of the interaction is 
screened by electrons of 3D gate electrodes. 

In order to calculate I-V curves we have derived boundary conditions at the contacts 
which take into account that the current leads of higher dimension play a role of heat bath 
for the quantum wire. 

Calculations of the transport properties are based on the equations of motion derived for 
the phase fields at the impurity site. The equations resemble equations of motion of a quan- 
tum pendulums. Therefore, there are solutions corresponding to rotation of the pendulum, 
that is solutions in which the phase increases linearly with time and the linear increase is 
accompanied by time-periodic oscillations. Since time derivative of the phase determines 
current in the wire, this corresponds to the regime of conductance in which dc current is 
accompanied by generation of ac current. The solution of the problem is complicated by 
strong fluctuations of the phase fields. For the case of spinless electrons we describe fluctu- 
ations under assumption that fluctuations are Gaussian. This assumption is justified in the 
limits of large voltages and/or strong interaction between electrons, otherwise it is a model 
assumption. In case of spinful interaction the problem becomes more complicated because 
in practically important cases K a is not small, therefore, fluctuations in the spin channel 
are neither small nor Gaussian. However, fluctuations in the spin channel can be treated 
strictly in case of spin-rotation invariant interaction, when interaction does not dependent 
on spin and K a = 1. Then we can use refermionization in the spin channel which makes the 
impurity Hamiltonian quadratic in spin variables. This enables us to solve the problem in 
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the limits of strong short-range interaction (K p <C 1) or moderate long-range interaction. 

As non-adiabatic contacts between the quantum wire and leads of higher dimensions pro- 
duce Friedel oscillations in the quantum wire, similar non-stationary effects in the electronic 
transport can be induced by non-ideal contacts. 

We have found the dynamical regime of conduction when the applied voltage exceeds the 
threshold value. The magnitude of the threshold voltage is determined by 2/c^-component of 
defect potential reduced by quantum fluctuations. Thermal fluctuations further reduce the 
threshold voltage, therefore, the effect of impurity on conductance is suppressed at tempera- 
tures T > To ~ Vt- Furthermore, fluctuations destroy the effect of the impurity in relatively 
short wires, shorter than the length of the order of v/Vt- Thus the effect considered in our 
work can be observed at low enough temperature in a relatively long quantum wire with a 
single impurity or a non-ideal contact, with minimal length and maximal temperature being 
related to magnitude of impurity potential and the strength of inter-electronic repulsion. 

The impurity potential W can be of different origin and of different strength. If the 
defect is induced by an impurity atom in the conduction channel then the potential can be 
quite large, of the order of the Fermi energy, but if the defect is made artificially, say, by a 
potential of a point contact, then W can be quite small. In semiconductor based quantum 
wires with shallow impurity W is expected to be of the order of few millivolts. In this case 
the range of frequencies of generated ac signal can be quite large, up to practically important 
teraherz region, depending on material of the quantum wire and the origin of the defect. 
We expect that one of the methods of experimental observation of current oscillations is to 
search for resonant Shapiro states in the presence of the applied ac-voltage. 

We have also found that I-V curves of the quantum wire with defect depend on the spin 
bias applied to the system. Thus spin bias and spin current can be measured electrically. 
This is related to violation of the spin-charge separation at the impurity site. In the consid- 
ered case of spin-rotation invariant interaction [K a = 1) the spin current cannot be induced 
by electric voltage if there is no spin bias. It is interesting to note, that in Gaussian model 
for spin channel one can obtain non-zero spin current as the result of spontaneous symmetry 
breaking, even in case when spin bias is absent. However, as the Gaussian model cannot 
be justified when K a is not small compared to 1, these results cannot be taken seriously. 
Nevertheless spin current can be induced in the absence of spin bias if the impurity is mag- 
netic^. In this case the electric voltage induces not only the electrical current, but the spin 
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current as well. 
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